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In the present time we observe a growing number of pub- 
lications where the, so-called, flow equations are successfully 
used to diagonalize Hamiltonians by means of an appropri- 
ate unitary transformation. Here we discuss and compare 
the flow equations (FE) method (proposed in 1994) with the 
method of one step continuous unitary transformations (OS 
CUT) (proposed in 1982). It is shown that the FE method 
can be considered as a generalization of the OS CUT approach 
to the case of parameter dependent generator. The OS CUT 
method gives linear differential equations for the diagonaliza- 
tion procedure. In the FE method the system of differential 
equations is nonlinear. Finally we discuss the generalization 
of idea of continuous unitary transformations for the case of 
quantum equations of motion (Heisenberg picture and density 
matrix) . 



I. INTRODUCTION 

Unitary transformations play an extremely important 
role in physics. For example, with the help of unitary 
transformation sometimes it is possible to simplify a 
problem of interacting quasiparticles by eliminating "in- 
convenient" interaction terms from initial Hamiltonian 
and construct corresponding effective interaction terms. 

In the present time we observe a growing number of 
publications! (seejl}, @, @, @, @, §, @, [§, §, @, 
0|, Q, @, Q J, where the, so-called, flow equations 
(FE) method is successfully used to diagonalize Hamilto- 
nians by means of an appropriate unitary transformation. 
This method of continuous unitary transformations looks 
very convenient and universal. The transformed Hamil- 
tonian appears as a result of solving differential equations 
on some formal parameter. 

It should be noted that quite the same idea and an 
analogous method (Safonov, 1982 @, 1983 @) to con- 
struct nonlinear unitary transformation was published 
far before the original paper on FE method (Wegner, 
1994 [|lj) and just did not attract enough attention of re- 
searches. For convenience, we shall call the approach of 
Refs. ]l5| ], 1 16 1 as a one-step continuous unitary transfor- 
mation (OS CUT) method. 

The goal of the present paper is to discuss and com- 
pare FE and OS CUT methods. As it will be shown 
later, the FE method is a generalization of OS CUT ap- 
proach. The only difference between the methods is that 



the unitary transformation generator is assumed to be 
dependent on a formal parameter in the FE method and 
the generator does not depend on the parameter in the 
OS CUT method. In the CUT method we obtain linear 
differential equations. In the FE method we should solve 
a nonlinear system of differential equations. Below we 
shall demonstrate both methods by several examples. 

It is interesting that the OS CUT method idea can be 
applied to the case of quantum equations of motion. Re- 
cently Mista and Filip |^3| have developed a method of 
non-perturbative solutions of nonlinear Heisenberg equa- 
tions. We shall briefly discuss this approach at the end 
of this paper. 



II. CONTINUOUS UNITARY 
TRANSFORMATION METHOD 

A. One step continuous unitary transformation 

Here we shall formulate the key idea of the OS CUT 
method 0I , |ll| . A general form of unitary transforma- 
tion can be written as 



H(0) 



He 



-8R 



(1) 



where is a formal parameter, i?is an anti-Hermitian op- 
erator = -R), H and H(6) are the initial and trans- 
formed Hamiltonians, respectively. The expression (|l]) is 
the solution of the equation 



>> 



[R, H(6)] 



(2) 
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with the initial condition H(Q) = H. In this form the 
Hamiltonian is "rotated in an operator space around i?." 

In order to solve this equation we should write the most 
general form of H{0) as the expansion in terms of opera- 
tor combinations with unknown ^-dependent coefficients. 
The most general form of R as anti-Hermitian operator 
combinations also should be used. 

After substituting these general (for H(9) and R) ex- 
pressions into (||), one can obtain a set of linear differ- 
ential equations by comparing coefficients in analogous 
operator compositions in both sides. Solving these equa- 
tions with the initial conditions, we obtain the trans- 
formed Hamiltonian H{0). In order to eliminate 'incon- 
venient' terms one needs to put their coefficients (for ex- 
ample, for = 1) equal to zero. This condition defines 
the choice of R. 

Usually it is simple to find the anti-Hermitian operator, 
which eliminates given non-diagonal term. The OS CUT 
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method first has been applied for the spin Hamiltonian 
diagonalization fl6|| . Then it was successfully used 
in physics of nuclear spin waves JITt , magnon-impurity 
interactions @, theory of superconductivity ||l9|] and for 
eliminating three-boson interactions p0[ . 



B. Flow equations for Hamiltonians 

The flow equations method |l| begins with a general 
form of a unitary transformation: 



H{£) = U{t)HU\l), 



(3) 



where £ is a continuous flow parameter. The Hamiltonian 
H{£) is transformed from the initial Hamiltonian H(0) = 
H. Differentiation (||) yields 



-H(t) = [r 1 (e),H(t)} 



with the generator 



r,(£) 



dU{£) 
d£ 



U\t) = -rf{£). 



(4) 



(5) 



The generator rj(£) should be chosen in such a way, that 
the off-diagonal matrix elements decay. A good choice is 
defined by 



r 1 (£) = [H d (£) ) H(£)} 



(6) 



where Hd{£) is the diagonal part of H{£) (comments FLQ). 

One can see that the FE method utilizes the same key 
tool (first proposed in the OS CUT method): the unitary 
transformation is represented as a solution of differential 
equations on a formal parameter. With an accuracy of 
notations the only difference between FE procedure and 
the OS CUT method is the parameter dependence of the 
generator -q{£). If -q does not depend on £ we have U(£) = 
exp(rj£). This case, obviously, exactly corresponds to the 
OS CUT method. If U{£) = exp(RF(£)), where R is 
the parameter independent anti-Hermitian operator and 
F{£) some function on £, then r)(£) = {dF(£)/d£)R. The 
Eq.(^) in this case can be written as 



(dF(£)/d£)d£ 



H{£) = [R,H{£)]. 



(7) 



This equation is reduced to Eq.(Q) using a simple 
change of variable dQ = (dF(£)/d£)d£. For U{£) = 
exp (RiFi(£)+R 2 F2(£)) with commuting i?i and R 2 anti- 
Hermitian operators one has n(£) — (dFi(£) / d£)R\ + 



1 A rigorous criterion that this choice is most optimal seems 
has not been obtained. 

2 The condition of diagonality and the condition that n{£) = 
[Hd{£), H(£)] — are not equivalent to each other. 



(dF 2 (£) / d£)R 2 . The diagonalization procedure in this 
case can be considered as a two transformations with 
U x {£) = exp(R 1 F 1 (£)) and U 2 {£) = exp (R 2 F 2 (£)) uni- 
tary operators. Each transformation is described by (Q) 
with R u Fx{£), or R 2 , F 2 (£). 

In the most general case the ^-dependent generator 
transformation can not be reduced to the parameter inde- 
pendent case. In such a situation the FE method can give 
principally different results than the OS CUT method. 
Let us consider several examples. 



III. QUADRATIC FORM DIAGONALIZATION 

Stein d} has considered the following quadratic Hamil- 
tonian: 



(8) 

where a q , a q and 6 q , 6 q are the creation and annihilation 
Bosc operators. Applying flow equation 



with 



-H 1 (£) = [r l0 (£),H 1 (£)] (9) 



%W = JE-9qW« 6 -q- a q^q)' ( 10 ) 



the following nonlinear differential equations have been 
obtained: 

| = | = -/q W3q (11) 



The solutions of (IllD are 



/ q W = e q coth(£e q + 4(q)), 



e q sgn(7 q ) 



sinh^e q + 4(q)) : 



(12) 



where e c 



(l- 7 2)V2 j£o(q)= i ln (i±£a) and 



e« = yj%(t)-9l(e) = yj%(0)-9l(0) (13) 

is the spectrum of the diagonal (at £ — > oc) Hamiltonian 
(|). 

Let us now consider the analogous diagonalization pro- 
cedure in the framework of the OS CUT method. For 
convenience, the same notations will be used. Simple 
analysis shows that the only anti-Hermitian operator that 
can eliminate a q £>! q + a q 6_ q in (^) is a q 6^L q — » q &- q - 
Thus, the generator (independent on £) should be taken 
in the form: 
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R = \ G q( a q 5 -q _ ft q 6 -q)' 



(14) 



where G q is unknown function independent on £. From 
dHi(£)/d£ — [R, Hi(£)] it is simple to obtain linear dif- 
ferential equations 



-/ q (£) = -G q5q W, - 5q (^ = -G q / q (£). (15) 
The solution of (|l^) is simple: 

m = / q (0)+ gq (0) e -^ + /q(0)- g q(0) eGq , ; (16) 
^ = /q(0)+gq(0) -g n < _ /q(0)-gq(0) cGo / 

Let the Hamiltonian (||) be diagonal at I = 1: <? q (l) = 0. 
From this condition we immediately find the unknown 
function G q : 

Gq ~2 ln U(0)- 5q (0)J 
and the spectrum 



eq = /q(l) = V/q(°)-5q(0) 



(18) 



Thus, we see that in the case of exact diagonalization of 
quadratic Hamiltonian the results of FE and OS CUT 
methods coincide with each other. 



IV. ELECTRON-PHONON COUPLING 

Consider now the elimination of electron-phonon inter- 
action to obtain the effective electron-electron scattering 
in the framework of OS CUT and FE methods. The 
Hamiltonian is (see, B): 



where 



H — Hq + H e - p h + H e _ e , 



(19) 



(20) 



describes free phonons (a q , a q ) and electrons (c k , c k ), 
H e _ ph {£) = M k , q WaL q 4 +qCk + h.c. (21) 

k,q 

describes the electron-phonon interaction and 

He-e(l) = £ Vk,k', q WCk +q 4_ q Ck'C k (22) 
k,k',q 

is the Hamiltonian of electron-electron interaction. The 
initial conditions arc: 



M k , q (* = 0) = M q> Vk,k',q(^ = 0)=Vk,k',q(Q). (23) 

Let us first consider the diagonalization (elimination 
of (pl|)) in the framework of OS CUT method. From a 
simple analysis it follows that the parameter independent 
generator should be taken in the form: 



R = £ ^k, q a- q c k + q c k - b.c. 

k,q 



(24) 



The diagonalization procedure gives the following linear 
differential equations: 



0M k , q (l) 
d£ 



-ak,qi?k,q 



(25) 



and for the electron pairs (with zero momentum) 

dVk,_ k ,qW 



d£ 



i?k,qM_k_q,q(£) - fi-k-q.q M k , q (I) . 

(26) 



From (|25| ) we obtain 

M k , q W = M k , q (0) - a k , q i?k, q l 



(27) 



Let the electron-phonon Hamiltonian ( |2l| ) vanishes at £ = 
1: M kiq (l) = 0. From this condition we immediately find 



-Rk.q = M kiq (0)/a kj q. 



(28) 



Simple solution of fl26|) with (|27|), (|28j ) and initial condi- 
tions (E3) gives 



U k ,-k,q(l) = Vk,_ k , q (0) 



W q - Ok+q - £k) 



(29) 



This result corresponds to the Frohlich's effective inter- 
action obtained in the second order of the perturbation 
theory (see, e.g., @). 

Let us now consider the FE method. The generator is 
taken in the form || : 

n{£) = [H Ql H^ ph (£)] (30) 
= £ Mk : q(f)ak, q a-qC k+q Ck - h.c. 

k,q 



with the energy difference 

«k,q = ^k+q — £k + W q . 



(31) 



This generator yields several contributions to dH(£) jdl = 
[r)(£),H(£)]. The contribution to the change of M k>q (£) 
results from [ry(^),ffo] 



Q£ - -a k qM kt q{£) 



with the solution 



M k , q (£) = M q exp(-a^ql), 



(32) 



(33) 
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where M q is the initial electron-phonon coupling. For 
c^k.q 7^ the solution (33) eliminates the electron-phonon 
interaction from the Hamiltonian at £ — > oo. The con- 
tribution to the effective electron-electron interaction is 
obtained from \r]{£), H e - P h(£)]. For the electron pairs 
(with zero momentum) we have the following nonlinear 
differential equation: 



d^k,-k,q(l) 

d£ 



-(a k , q + a_ k _ q , q )M k>q (£)M_ k _ q , q (£) 

(34) 



with the solution 



k,q 



(oo) = Vk,-k, q (0) - 



M^ q 



(£k+q - £k) 5 



(35) 



This result differs from (|29| ) (plus sign between the two 
squares in the denominator). The difference is explained 
as: "perturbation theory for Hamiltonians is not uniquely 
defined" §. 

In the FE method a choice of r){£) is motivated to be 
defined by [Ho, H e - P h(£)] (|3(]). However, a rigorous cri- 
terion that this choice (or, another one) is most optimal 
for a given problem seems has not yet been formulated. 



V. THREE-BOSON HAMILTONIAN 

In this example we consider the elimination of forbid- 
den three-boson terms from the Hamiltonian 



H(£) = ^£k&]A + H 3 (l) + H 4 (£), 



(36) 



where 



k,q 



and 



Ha(£) = \ $ k 1 ,k 2 ;k 3 ,k 4 W6 Ul 6 U2 6k 3 6k 4 

ki,k 2 ;k 3 ,k 4 

xA(ki +k 2 -k 3 -k 4 ). 



(38) 



Here £ is a formal parameter and A(k) is the Kronecker 
delta function. For simplicity, we shall neglect the energy 
£k dependence on £. 

Using OS CUT method, we write the parameter- 
independent operator in the form |20|: 



k,q 



(39) 



Calculating the commutators [R, ^ k £k&k^k] and 
[R, H%{£)], we obtain the following linear differential 
equations: 



^*k,qW = -3(£k + £q + £-k-q)i?k,q, (40) 



df 



^k 1 ,k 2; k3,k 4 W = -6i?k 1 ,k 2 * k3>k4 W 



(41) 



Solving these equations with the initial conditions 
*k!,k 2 (^ = 0) = * kl ,k 2 , ^k u k 2 -MM^ = 0) = 0, and 
demanding v J , ki,k 2 (^ = 1) = 0, we obtain 



i?k.c 



k.q 



3(£k 



£-k- 



(42) 



and 



$ kl ,k 2 ;k 3 ,k 4 (l) = -^.k.*^ 

1 



1 



+ - 



£ki + £k 2 + £ -ki-k 2 

(43) 



£k 3 + e k 4 + £-k 3 -k 4 

It should be noted that this result for the effective four- 
boson interaction coincides (with the additional condi- 
tion = with the result of classical diagonaliza- 
tion of ([37|) (when and 6k are c-numbers) by canonical 
quasi-linear transformation j^TJ]. 

Let us now consider the same problem in the frame- 
work of FE method. The (parameter-dependent) gener- 
ator in this case is: 



VV) = I £{A<, q *k, q W&k&q& t -k- q - h.c.}, 



k,q 



where 



/?k,q = £k + £q + £-k-q- 



(44) 



(45) 



Corresponding differential equation for the three-boson 
amplitude has the form: 

^*k,qW = W) (46) 

with the solution (compare with ([33])): 

*k,qW = *k, q exp(-/3^). (47) 

The differential equation for the four-boson effective 
amplitude in this case is 



rf$kl,k 2 ;k 3 ,k 4 (l) 

d£ 



-2(/3k 1 ,k 2 +/3k 3 ,k 4 )*k 1 ,k 2 W* k3 ,k 4 W- 

(48) 

The solution at £ — > 00 (when \& k , q (^) vanishes) has the 
form 



ki,k 2 ;k 3 ,k, 



(1) 



-2* kl ,k 2 * k3 ,k 4 ^2 



Aci,k 2 + Ac 3 ,k 4 



(49) 



B 2 

f - ki,k 2 ■ <-k 3 ,k 4 

We see that this formula differs from the OS CUT re- 
sult ( f43| ) and therefore ( pTjf ) does not correspond to the 
classical case 1231. 



4 



VI. NON-PERTURBATIVE SOLUTION OF 
QUANTUM EQUATIONS 

A. Heisenberg picture 

The Heisenberg equation for an operator A(t) 



This example can demonstrate how the order stored in 
magnetization is transformed to higher-order spin corre- 
lations. 

The author wish to thank G. S. Uhrig and A. K. 
Khitrin for valuable comments. 



ih-A(t) = [A(t),H(t)} 



(50) 



can be rewritten in the form: 



-A(t) = [iH(t)/h,A(t)}. 



(51) 



Formally this equation is similar to (Q) (or, (Q) for time- 
dependent Hamiltonian): time t plays a role of a formal 
parameter and iH(t)/h is an anti-Hcrmitian operator. 
As in the CUT method, we can write the most general 
form of A(t) as the expansion in terms of operator combi- 
nations with unknown time-dependent coefficients. After 
substituting these form into ( pl| ) with a given iH(t)/h, 
one can obtain a set of differential equations by com- 
paring coefficients in analogous operator compositions 
in both sides. These equations will be linear for time- 
independent H and nonlinear for time-dependent H(t). 
Then we solve these equations with the initial conditions 
defined by A(t) = A(0). Recently Mista and Filip §| 
have developed this idea as a method of non-perturbative 
solutions of nonlinear Heisenberg equations with an illus- 
trative example of two coupled harmonic oscillators. 



B. Density matrix 

It is obvious that a quite analogous idea can be applied 
for the density matrix. In this case we can write the 
density matrix equation as 



dt 



p(t) = [-iH(t)/h,p(t)}, 



(52) 



where —iH(t)/h plays a role of "transformation genera- 
tor". 

As an illustration consider a system of nuclear spins 
lj with a Hamiltonian H = Hq + Hi nt , where Hq = 
pB YU Ij is the Zeeman energy and H int describes spin- 
spin interactions. As an initial state assume the trans- 
verse magnetization. In the high temperature approx- 
imation one can write p(Q) = 1 — aj^j^j- Then 
the density matrix can be represented as p(t) = 1 — 
^^ =1 A/9„(t), where Ap n (t) describes compositions of 
n spins with time dependent coefficients, e.g., Ap±(t) = 
Ej i x j{t)Ij + Yj(t)Ij)- From Wo, -Api(t)] we obtain 
differential equations dXj (t) /dt = —tooYj (t) , dYj (t) /dt = 
u>oXj(t) and solutions Xj(t) — acosuiot, Yj(t) — 
asinajpt, where loq — pB/h is the nuclear magnetic reso- 
nance frequency Contributions to two spin motions can 
be found from [Ho, —Ap 2 (t)] + [H in t, — Api(t)] and so on. 
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